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SIGNAL  COHERENCE  RECOVERY  USING  ACOUSTO-OPTIC 
FOURIER  TRANSFORM  ARCHITECTURES 


INTRODUCTION 

The  program  in  acousto-optics  (AO)  conducted  by  the  X-ray  Astronomy  Branch,  Space  Science 
Division,  at  the  Naval  Research  Laboratory  (NRL)  has  a  special  astronomical  orientation- 
development  of  ultrafast  analog  signal  processors  for  astronomical  purposes.  Time-series  analysis 
provided  the  initial  development  challenges.  This  has  meant  that  special  emphasis  was  placed  on  sys¬ 
tems  that  would  analyze  periodic  and  nearly  periodic  functions.  Processers  of  this  type  also  have 
nonastronomical  applications.  Our  particular  approach  to  time-series  analysis  is  called  coherence 
recovery  (CR).  It  applies  generally  to  analyses  of  faint,  frequency-modulated  functions.  The  fre¬ 
quency  modulation  is  removed  by  a  trial-and-error  process  that  involves  extensive  computation.  The 
faint  signal  becomes  detectable  and  the  characteristics  of  the  frequency  modulation  are  simultaneously 
established. 

Overview 

In  last  decade,  the  maturing  technologies  of  AO  components  have  led  to  declining  costs,  making 
an  aggressive  development  program  for  fast  processors  economically  feasible.  A  two-dimensional  (2- 
D)  Fourier  processor,  Mark  II,  for  which  design  and  procurement  stages  have  been  successfully  com¬ 
pleted,  should  have  the  equivalent  throughput  of  two  Cray-2s  on  algorithms  such  as  very  long  Fourier 
transforms.  The  1989  effort  consisted  of  building  Mark  II  and  researching  adaptations  of  the  2-D 
design  for  digital  optical  processing.  Mark  ITs  key  assets  (two  Cray-2s,  inexpensive  to  use,  special 
features  to  isolate  chirped  signals)  make  it  attractive  for  a  range  of  applications. 

Motivation 

The  Space  Science  Division  AO  effort  is  unique  among  astronomical  groups,  and  the  astronomi¬ 
cal  application  is  unique  among  optical  computing  groups.  The  primary  thrust  of  the  program  has 
been  to  seek  improvements  in  techniques  involving  Fourier  transforms.  The  Mark  II  design  is  at  the 
frontier  of  the  current  development.  A  similar  processor  was  outlined  in  the  literature,  but  it  was 
never  constructed.  The  current  design,  however,  matches  recently  matured  AO  technology  with  real¬ 
time  astronomical  data  processing  demands.  A  description  of  one  astronomical  problem  for  which 
Mark  II  provides  a  solution  will  illuminate  analogous  applications  in  other  fields. 

This  has  meant  that  special  emphasis  was  placed  on  systems  that  would  analyze  periodic  and 
nearly  periodic  functions.  Processers  of  this  type  also  have  nonastronomical  applications.  Our  partic¬ 
ular  approach  to  time-series  analysis  is  called  coherence  recovery  (CR).  It  applies  generally  to  ana¬ 
lyses  of  faint,  frequency-modulated  functions.  The  frequency  modulation  is  removed  by  a  trial-and- 
error  process  that  involves  extensive  computation.  The  faint  signal  becomes  detectable  and  the 
characteristics  of  the  frequency  modulation  are  simultaneously  established. 

Manuscript  approved  April  8.  1989. 
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One  of  our  original  astronomical  motivations  was  the  search  for  millisecond  pulsations  from 
binary  stellar  systems.  Such  millisecond  pulsars  would  have  significant  ramifications  for  gravitational 
physics.  Reference  1  provides  a  scenario  whereby  a  neutron  star  can  be  spun  up  until  it  becomes  sub¬ 
ject  to  gravitational  radiation  reaction  instabilities.  The  object  then  becomes  simultaneously  a  pulsar 
in  X  rays  and  gravitational  waves.  The  envisioned  detection  strategy  is  to  find  the  signal  in  X  rays 
and  then  build  a  special  gravity  wave  antenna  with  narrow-band  response  tuned  to  the  frequency 
determined  in  the  X-ray  band. 

Because  the  pulsations  are  expected  to  be  weak  (<0.1%  modulation),  long  transforms  are 
required  to  achieve  sensitivity.  In  addition,  the  detection  of  such  pulsations,  in  both  electromagnetic 
and  gravitational  radiations,  requires  that  the  effect  of  orbital  motion  be  removed  before  performing  a 
frequency  search  by  means  of  a  Fourier  transform  (in  the  binary  systems,  which  are  likely  millisecond 
pulsar  candidates,  the  orbital  parameters  are  partially  or  totally  unknown  but  constrained).  Otherwise, 
loss  of  coherence  (i.e.,  frequency  smearing)  renders  a  weak  signal  unobservable.  The  approximate 
removal  of  orbital  motion,  termed  coherence  recovery  (CR)  (Fig.  1),  can  be  effected  by  performing  a 
one-parameter  grid  of  quadratric  time  transformations  and  subsequent  Fourier  transforms  [2,3];  one 
transformation  of  the  set  will  provide  optimal  recovery  over  the  necessary  long  integration  time.  The 
ideal  utilization  of  CR  is  in  conjunction  with  a  very  large  collecting  aperture.  In  Ref.  4  we  proposed 
the  construction  of  a  100-m2  array  that  will  detect  pulsations  with  amplitudes  <0.01%  (Fig.  2). 

The  computational  burden  required  for  CR  is  determined  by  the  large  number  of  trial  time 
transformations  (proportional  to  the  square  of  integration  time)  required  to  effect  recovery.  Digital 
implementation  in  real  time,  required  for  gravitational  wave  searches  and  desirable  for  some  onboard 
spacecraft  data  processing,  can  be  prohibitively  expensive,  requiring  as  many  as  10  to  100  Cray- 
equivalents  of  computational  resource.  Acousto-optic  Fourier  processors  provide  direct  analog  tech¬ 
niques  for  fast,  real-time,  second-order  time  transformations  and,  hence,  the  best  possible  CR.  Thus 
they  are  well-matched  to  various  problems  that  involve  rapid  searches  for  frequency-smeared  signals. 


Ppulse:  20  ms  Porb:  3800  s  A:  2.35  1-s  Ntrans:  39 
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Fig.  I  —  Digital  simulation  of  coherence  recovery  (a)  coherent  signal  modulated  sinusoidally,  plus  noise  and 
(b)  optimum  trial  transformation  recovers  -90%  of  original  signal  strength. 


2 


NRL  REPORT  9214 


HfcAOA-1:  1m2  (1977) 


SINGLE  ARRAY  MODULE 


TRUNNION-MOUNT 
POINTING  CONTROL 


XLA  ARRAY  (64  MODULES.  100  m2) 


Fig.  2  —  The  X-ray  large  array 


The  main  advantages  of  acousto-optics  are  high  throughput  at  low  cost  and  the  promise  of  appli¬ 
cation  in  orbit  to  tasks  such  as  onboard  timing  analysis.  On  board  applications  have  the  additional  vir¬ 
tue  of  being  resistant  to  soft  upsets  caused  by  charged  particles.  Steady  improvement  in  ancillary 
technologies,  such  as  chirp  generators,  is  another  attraction;  this  means  that  ambitious  processors  can 
now  be  developed  at  low  cost. 

Two-dimensional  acousto-optic  processors  have  many  uses — not  only  on  orbital  platforms  but 
also  in  ground-based  applications.  In  addition  to  fast  Fourier  transforms  (FFTs),  acousto-optic  imple¬ 
mentations  of  Hilbert  transforms,  ambiguity  functions,  general  correlation,  and  convolution  are 
directly  realized  with  analog  techniques,  thereby  greatly  surpassing  the  speed  of  traditional  electronic 
devices.  Development  of  acousto-optic  processors  for  “speed  of  light”  matrix  manipulations  using 
n-ary  coding  is  rapidly  progressing.  Many  astronomical  data  processing  applications  exist  for  these 
near-real-time  digital  optical  processors.  The  Space  Science  Division’s  acousto-optic  program,  with 
emphasis  on  astronomical  applications,  places  NRL  in  a  unique  role  in  both  the  astronomical  and  the 
optical  computing  communities. 
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(b)  2-D  Fourier  processor  under  construction 
Fig.  3  —  NRL  acousto-optic  computer  project 
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History  and  Accomplishments 

The  program  has  been  active  since  1986;  during  this  time  the  Mark  I  system  was  constructed 
and  tested,  and  Mark  II  was  designed.  This  is  summarized  in  Fig.  3.  Mark  I  is  a  time-integrating 
architecture  of  modest  capability  (Fig.  4).  It  is  one-dimensional,  in  the  sense  that  the  readout  format 
is  a  1-D  display  of  intensity  (representing  Fourier  power)  vs  position  (representing  frequency).  After 
initial  demonstration  of  system  concept,  Mark  I  was  upgraded  to  increase  its  speed  by  a  factor  of  10 
and  to  increase  the  number  of  frequency  channels  (Fig.  5).  The  scheme  for  CR,  described  briefly 
above,  which  had  earlier  been  developed  in  software  for  pulsar-searching,  was  found  to  have  a 
hardware  realization  that  could  be  achieved  with  the  addition  of  a  tunable  chirp.  This  technique  is 
general,  and  it  carries  over  to  Mark  II.  A  paper  describing  the  CR  technique  and  its  acousto-optical 
implementation  was  presented  at  a  meeting  of  the  American  Astronomical  Society  [5].  A  paper  on 
the  same  subject  was  presented  at  a  symposium  of  the  International  Society  of  Optical  Engineering 
[6]. 


Mark  II  design  also  occupied  much  of  1987.  The  goal  is  a  very  fast  processor  that  can  provide  a 
2-D  output  format  having  greatly  increased  bandwidth  and  throughput.  A  design  was  selected  that 
can  be  realized  in  various  forms,  each  suitable  for  a  specific  application.  We  began  with  a  simple 
version  involving  only  limited  additional  purchases  and  considerable  recycling  of  Mark  I  components. 
This  design  was  verified  by  using  symbolic  manipulator  programs  to  check  the  integral  transforms 
representing  the  nominal  performance.  A  key  component  of  Mark  II,  the  fast  chirp  driver  that  imple¬ 
ments  coarse  separation  of  frequencies,  represents  an  important  improvement  in  2-D  Fourier  proces¬ 
sors. 


KNIFE 

FILTER  BEAM 


Fig.  4  —  One-dimensional  time-integrating  Fourier  transform  architecture 


5 


MARK  I  ARCHITECTURE:  ONE-DIM  FIRST  FOURIER  TRANSFORMS  PRODUCED  BY  OPTIMIZATION  IMPROVES  RESOLUTION, 
T-I  FOURIER  TRANSFORM  PROCESSOR  MARK  I.  4  KHZ  AND  8  KHZ  SIGNALS.  SPEED,  DYNAMIC  RANGE.  10  KHZ  AND 
BUILT  APR  -  MAY  ’88  NO  OPTIMIZATION.  JULY  *88  IS  KHZ  SIGNALS.  SEP  -  OCT  ’86. 


WOOD.  NORRIS,  AND  SMATHERS,  JR 


El  3TC 

ES  <"X 

>} 

T  i 

>1  If  . 

si 

"  f  N.  • 

§  Sr 

~  Mr 

>  \\ 

o  U 

\\  1 
l  iv 

CO  rt 

i  \i 
>'*• 

it  * 

I  J. 

I 

_ _ _  .  _  ,i  r _ 

UJ 

<r 

O 

• 

< 

z 

O 

3 

< 

M 

O' 

K 

CO 

-J 

< 

< 

X 

to 

H 

UJ 

Z 

K 

to 

UJ 

U 

o 

2 

UJ 

Q- 

< 

a: 

2 

a 

oc 

O 

z 

o 

o 

3 

u 

UJ 

U. 

o 

z 

z 

M 

• 

HI 

to 

O 

a: 

HI 

• 

UJ 

z 

a. 

* 

o 

o 

2 

a: 

Q- 

OL 

UJ 

< 

M 

•  » 

X 

cc 

UJ 

D 

> 

o 

o 

< 

(r. 

LL 

* 

CO 

8  8 
CO  t- 

H3MOd  aaznvwyoN 


UJ 

O' 

CM 

• 

D 

® 

o 

00 

W 

u. 

» 

z 

> 

CO 

o 

o 

M 

z 

o 

h- 

M 

> 

_J 

• 

UJ 

o 

z 

a: 

CO 

< 

a. 

UJ 

X 

a: 

o 

Ur 

o 

UJ 

>- 

»- 

u 

2 

< 

z 

D 

2 

UJ 

X 

M 

3 

H 

K 

O' 

O 

-1 

UJ 

UJ 

D 

a: 

a. 

u. 

<0 

UJ 

N. 

> 

to 

oc 

UJ 

-J 

UJ 

M 

UJ 

X 

o 

U 

M 

a. 

< 

a. 

O 

to 

cc 

UJ 

a 

-J 

X 

• 

CO 

u 

CO 

< 

< 

so 

z 

CD 

* 

UJ 

1 

> 

oc 

o 

UJ 

z 

H 

-J 

o 

< 

a. 

z 

• 

2 

o 

2 

O 

M 

D 

u 

to 

a: 

►— 

K 

u 

< 

• 

UJ 

o 

a. 

2 

z 

to 

00 

M 

M 

U 

UJ 

z 

a 

O 

UJ 

D 

K 

o: 

t- 

UJ 

M 

* 

u- 

-J 

Z 

u. 

0l 

HI 

M 

2 

-J 

a 

< 

8  8  8 

d3MOd  CJ3ZnvWdON 


8 


"  x 

CD  ~ 
^  > 
o 

uo  Z 
uj 

3 

CM  O 
LU 
DC 


O  ^  OO  CM  (D 

CO  CM  r“ 

aanindwv  aiznvwuoN 


(A 


& 


'£ 


6 


ACOUSTO-OPTIC  FOURIER  TRANSFORM  PROCESSORS 


NRL  REPORT  9214 


MARK  I:  ONE-DIMENSIONAL  TIME-INTEGRATING  FOURIER  PROCESSOR 
Design  and  Concept 

The  design  of  Mark  I  resurrects  a  proven  architecture  for  exploratory  purposes.  The  original 
version  of  this  time-integrating  Fourier  architecture  was  first  proposed  and  implemented  by  the  Opti¬ 
cal  Sciences  Division  at  NRL  [7],  The  use  of  a  similar  architecture  provided  a  working  model  to 
implement  major  modifications  (CR)  and  advance  to  a  design  for  a  2-D  Fourier  system 

Figure  4  shows  the  geometry  of  the  Mark  I  processor;  this  may  facilitate  comprehension  of  the 
operating  principles  involved.  Mark  I  is  composed  of  a  front  end,  where  data  are  applied  to  an  AO 
modulator  to  modulate  a  laser;  a  cylindrical  telescope  operated  in  reverse  that  produces  a  horizontal 
sheet  beam;  a  central  interferometer  segment  that  creates  a  Fourier  kernel;  and  a  time-integrating  out¬ 
put  arm  in  which  the  kernel  beats  with  the  (temporal  domain)  input  to  compute  the  (frequency 
domain)  output.  The  kernel  is  actually  realized  upon  recombination  of  the  two  beams  at  the  detector 
camera  plane,  where  most  extraneous  phase  terms  cancel  identically  Here,  the  spatial  coordinate  of 
a  1-D  sensor  is  related  to  the  Fourier  frequency  coordinate.  Stops  designed  to  remove  the  (unwanted) 
undiffracted  parts  of  the  beam  also  appear.  Reference  8  provides  a  comprehensive  treatment  of  many 
aspects  of  AO  signal  processing. 

In  AO  Fourier  architectures,  the  signal  may  be  introduced  by  intensity  or  frequency  modulation 
of  the  light  beam;  for  the  design  shown  in  Fig.  4,  the  method  is  intensity  modulation.  The  signal  is 
subsequently  carried  in  both  legs  of  the  interferometer.  Here  two  counterpropagating,  chirped 
sinusoids  (of  form  exp  |/ar2),  a  =  constant)  are  impressed  onto  the  light  beams  by  frequency  modu¬ 
lation  using  Bragg  cells  (acousto-optic  beam  deflectors,  AOBDs).  The  momentum-matching  conditions 
between  the  light  and  sound  waves  within  the  AOBD  require  either  increasing  or  decreasing  the  light 
frequency  by  the  sound  frequency;  this  is  sometimes  referred  to  as  up  and  down  (frequency)  shifting 
of  the  light.  Selection  of  either  the  first  plus  or  minus  diffraction  order  determines  whether  the  light 
is  up-  or  down-shifted,  respectively.  The  Fourier  kernel  results  from  heterodyning  the  up-  and 
down-shifted  beams.  The  relative  light  frequency  lag  between  the  two  sheet  beams  (a  function  of  AO 
ceil  spatial  coordinate,  and  therefore  of  detector  coordinate)  provides  a  difference  frequency— the 
Fourier  kernel. 

It  may  help  to  think  of  this  arrangement  as  a  bank  of  frequency  filters.  One  such  filter  tests  for 
the  presence  of  a  periodic  component  at  or.e  discrete  frequency,  at  one  spatial  location  in  the  detector 
plane,  by  maintaining  a  (temporal)  sinusoid  that  is  heterodyned  with  any  arbitrary  input.  The  bank  of 
filters  maintains  a  series  of  such  sine  waves,  one  at  every  frequency  within  a  certain  bandwidth,  and 
the  frequency  so  maintained  advances  uniformly  as  a  function  of  spatial  position  within  the  device. 
At  the  output  camera,  the  detected  signal  builds  up  continuously  in  one  channel  over  the  course  of  an 
integration  cycle.  These  are  the  essential  features  of  a  time- integrating  Fourier  architecture.  In  addi¬ 
tion,  heterodyning  of  the  input  with  all  the  filters  whose  frequencies  differ  from  the  input  signal  fre¬ 
quency  produces  a  light  bias  that  builds  up  in  all  channels;  the  signal  is  detected  against  this  back¬ 
ground. 

By  contrast,  a  space-integrating  Fourier  architecture  simply  uses  one  Bragg  cell  to  diffract  an 
input  beam  at  angles  that  are  proportional  to  the  frequencies  present  in  the  signal.  The  diffracted 
beam  is  then  passed  through  a  spherical  lens  that  accomplishes  the  Fourier  transform.  The  light  bias 
built  up  over  an  integration  period  is  insignificant  compared  to  that  in  a  time- integrating  processor  and 
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is  essentially  attributable  to  stray  laser  light.  Hence,  the  detected  light  represents  only  the  signal;  this 
is  an  important  feature  of  space-integration.  However,  space-integration  does  not  allow  for  extended 
temporal  integration.  Also,  in  the  context  of  AO  Fourier  processing,  a  2-D,  purely  space-integrating 
architecture  cannot  be  implemented.  The  distinctions  between  time-integrating  and  space-integrating 
architectures  must  be  kept  in  mind  for  the  discussion  of  Mark  n,  because  both  time-  and  space- 
integrating  geometry  are  combined  to  generate  the  2-D  folded  Fourier  output  of  that  architecture, 
while  the  beneficial  aspects  of  the  two  approaches  are  retained. 

Optical  Bench  and  Associated  Electronics 

The  project  moved  very  quickly  from  inception  to  development  of  Mark  I.  This  was  partly 
because  components  were  initially  borrowed  from  another  NRL  code.  Availability  of  those  com¬ 
ponents  dictated  details  of  the  first-generation  Mark  I  layout,  but  upgrades  were  later  made  to 
improve  performance. 

Upgrades  in  AO  equipment  during  the  Mark  I  phase  were  made  to  increase  the  data-processing 
speed  of  the  system.  A  new  AO  modulator  was  purchased  that  introduced  data  by  intensity  modula¬ 
tion  at  rates  as  high  as  10  MHz.  To  Fourier  analyze  data  introduced  at  these  higher  rates,  a  faster 
chirp  generator  was  obtained.  This  chirper  is  a  linear  voltage-controlled  oscillator  (VCO)  that  we 
were  able  to  operate  at  a  1-ms  chirp  rate  (the  effective  integration  period).  Faster  usable  chirp  rates, 
which  would  take  advantage  of  the  ultimate  speed  of  the  AO  modulator,  are  obtainable  from  the  VCO 
chirper  by  introducing  more  sophisticated  synchronizing  electronics.  More  important  from  the  stand¬ 
point  of  new  algorithms,  by  using  both  our  analog  chirper  and  a  digital  chirp  synthesizer,  we  were 
able  to  demonstrate  CR  in  the  AO  system  (described  below). 

An  IBM  PC  was  set  up  to  direct  control  of  the  AO  processor’s  input  and  output.  Computer  pro¬ 
grams  were  devised  that  communicated  with  a  digitizing  oscilloscope,  a  chirp  synthesizer,  and  a  data 
generator;  they  also  processed  the  signal  received  from  the  AO  system.  At  the  later  stages,  the  com¬ 
puter  provided  test  signal  inputs  for  the  data  generator.  These  inputs  were  tailored  to  simulate 
frequency-smeared  signals  on  which  the  AO  system  performed  simultaneous  CR  and  Fourier  analysis. 
The  next  two  sections  describe  the  progress  made  with  component  upgrades  and  the  invention  of  an 
AO  system  capable  of  CR. 

Performance 

Mark  I  underwent  a  series  of  improvements  from  July  1986  to  January  1987,  reflecting  our 
learning  curve  and  increasing  skills  in  applying  the  basic  principles  of  optical  techniques.  The  first 
test  of  Mark  I  was  very  crude;  initially  it  was  possible  to  resolve  an  8  kHz  input  band  into  a  few  fre¬ 
quency  channels.  Various  optimizations  improved  the  figures  of  merit  until  the  appearance  of  the 
analog  signal  was  a  sharp  spike.  Following  this,  much  further  improvement  was  realized  by  utilizing 
the  postprocessing  capability  of  the  PC  for  such  tasks  as  digitizing  output  signal  waveforms,  comput¬ 
ing  signal  minus  background,  and  finally  presenting  power  spectra.  By  December  1986,  it  was  possi¬ 
ble  to  exhibit  a  square  wave  input  analyzed  into  the  first  and  third  harmonics  in  the  proper  ratio. 
Design  iterations,  which  incorporated  faster  analog  input  electronics,  led  to  system  performance  closer 
to  the  intrinsic  limits  of  the  AO  components.  During  this  period  the  resolution  increased  from  <50 
channels  to  -400  channels  while  the  input  rate  increased  from  15  to  350  kHz.  Progress  in  the 
development  of  Mark  I  is  illustrated  by  the  quality  of  spectra  obtained  (Fig.  5). 
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Coherence  Recovery 

A  by-product  of  the  Mark  I  effort  was  the  discovery  of  an  AO  method  for  handling  sinusoids 
chirped  by  arbitrary  functions.  It  was  motivated  by  the  original  astronomical  application— to  detect 
faint,  fast  pulsars  in  binary  systems.  The  problem  is  that  the  binary  motion  of  the  source  alternately 
blue  shifts  and  red  shifts  the  signal,  i.  e.,  the  signal  is  frequency  modulated  (smeared  in  frequency 
space).  The  simplest  way  to  undo  this  process  would  be  to  compensate  the  data  by  using  the  binary 
orbit  parameters,  which  is  the  standard  technique  in  pulsar  work.  Unfortunately,  the  orbit  is 
unknown  in  nearly  all  cases.  A  full  search  of  the  space  spanned  by  the  amplitude,  frequency,  and 
phase  of  the  unknown  orbit  is  computationally  prohibitive  with  present  resources.  However,  searches 
of  a  one-parameter  grid  of  parabolae  (quadratic  transformations)  have  been  used  to  achieve  compara¬ 
tively  high  sensitivity  by  using  the  SSD  VAX  [2].  The  set  of  parabolae  can  be  thought  of  as  a  grid  of 
approximations  to  the  actual,  unknown  orbit.  (The  VAX  takes  —1/4  h  to  do  a  220-point  Fourier 
transform.  An  optimized  Connection  Machine  FFi  will  take  —  2  s  for  220  points  [9],  This  is  part  of 
the  overall  program  motivation— AO  processors  would  be  much  faster.)  This  algorithm  may  be  used 
to  achieve  approximate  CR  over  an  interval  that  is  short  compared  to  the  orbital  period.  Our  simula¬ 
tions  (Fig.  1)  show  that  >90%  of  the  power  actually  present  in  the  coherent  signal  can  be  recovered. 

We  have  discovered  an  analog  AO  technique  for  CR  that  is  identical  to  a  quadratic  transforma¬ 
tion.  The  technique  was  implemented  with  minor  additional  hardware  development  [6],  We  found 
that  CR  can  be  achieved  in  Mark  I  by  making  the  two  chirps  slightly  unequal.  Varying  the  difference 
in  chirp  accelerations  over  a  small  range  generates  the  required  set  of  quadratic  transformations.  The 
resulting  modification  to  the  Fourier  kernel  includes  an  additional  fringe  term  (with  entirely  negligible 
effect  in  practical  applications)  and  a  quadratic  phase  term  that  accomplishes  the  exact  function  of  the 
grid  of  parabolae  as  in  the  software  implementation.  Although  many  conceivable  applications  involve 
sinusoidal  (or  arbitrary  functional)  modulation  over  a  fraction  of  a  cycle,  for  which  approximate  CR 
is  obtainable,  linear  frequency  excursions  (red  or  blue,  growing  linearly  with  time)  can  be  -  100% 
recovered  with  this  algorithm.  Linear  frequency  excursions  >25%  have  been  successfully  recovered 
in  Mark  I.  The  fact  that  the  technique  is  analog  (no  complex  digital  premodification  of  chirp 
accelerations  is  necessary)  makes  it  fast  and  very  attractive.  Certain  specific  applications  might  more 
profitably  search  a  different  family  of  functions  to  effect  a  good  coherence  recovery.  In  this  case, 
preparation  of  chirps  might  be  implemented  with  dedicated  digital  hardware,  rather  than  with 
software.  The  quadratic  transformation  often  will  be  the  best  one-parameter  approximation.  Appendix 
A  develops  mathematics  for  CR  for  the  case  of  pulsars  in  binary  orbits,  together  with  AO  realization 
of  this  case. 

Astronomical  applications  of  the  CR  algorithm  include  real-time  searches  for  radio  and  X-ray 
pulsars,  searches  for  continuous-wave  gravitational  radiation  from  the  most  rapidly  spinning  pulsars, 
and  SETI  (search  for  extraterrestrial  intelligence).  In  addition  to  the  astronomical  uses,  potential 
applications  of  the  CR  algorithm  implemented  with  AO  processors  exist  in  other  fields,  i.e.,  in  prob¬ 
lems  where  the  sensitivity  for  real-time  detection  of  chirped  signals  may  be  increased  by  compensat¬ 
ing  the  effect  of  modulation  by  an  arbitrary  function.  The  algorithm  is  applicable  to  situations  where, 
during  the  observation  interval,  the  phase  (frequency)  modulation  is  sufficiently  close  to  quadratic 
(linear)  in  time,  or  where  the  observation  interval  is  limited  so  that  this  condition  obtains.  Such  dev¬ 
ices  may  be  adapted  to  problems  involving  propagation  of  waves  in  inhomogeneous  media,  e.  g.,  in 
such  disciplines  as  seismology  and  underwater  acoustic  signal  reception. 

With  the  emergence  of  the  CR  concept,  Mark  I  development  was  regarded  as  mature.  A  paper 
was  published  describing  the  general  CR  concept  and  its  realization  in  a  1-D  time-integrating  design 
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[6],  At  this  point  in  the  program,  limitations  were  being  approached  because  of  the  1-D  nature  of  the 
Mark  I  device.  It  would  be  possible  to  improve  upon  the  Mark  I  and  make  small  packages  suitable 
for  various  purposes.  However,  in  such  systems  integration  time  (and  therefore  sensitivity)  is  limited. 
Also,  the  1-D  format  affords  intrinsically  low  resolution  and  therefore  comparatively  poor  perfor¬ 
mance  when  measured  against  a  computer  such  as  a  Cray.  The  major  purpose  of  the  program  was  to 
build  simple  systems  that  could  compete  with  Crays  on  algorithms  such  as  Fourier  transforms. 

MARK  II:  HYBRID  (SPACE  AND  TIME  INTEGRATING)  TWO-DIMENSIONAL  FOURIER 
PROCESSOR 

Mark  II  is  a  2-D  hybrid  architecture  with  performance  figures  of  merit  that,  in  some  respects, 
exceed  those  of  a  Cray -2.  The  2-D  output  format  resembles  a  raster  scan  of  frequency  space. 
Hybrid  refers  to  the  combination  of  space- integrating  and  time-integrating  subsystems  that  are  used  io 
achieve  an  effective  one-parameter  (frequency)  readout  format.  The  entire  Fourier  transform  is  real¬ 
ized  in  piece,  on  repeated  passes  of  the  input,  more  rapidly  than  a  Cray  can  accomplish  the  same 
computation.  The  throughputs  (bits/second  processed  and  recorded)  are  comparable.  The  AO 
processor’s  bit  accuracy  is  less  than  that  of  the  digital  computation,  so  that  the  principal  gain  is  in 
speed.  In  fact,  Mark  II  computes  a  106-point  Fourier  transform  approximately  500  times  faster  than 
does  a  Cray -2. 

A  significant  compromise  has  been  accepted  in  the  area  of  greatest  advantage,  namely  speed. 
This  compromise  is  the  repeated  introduction  of  the  input,  made  necessary  by  the  choice  of  a  com¬ 
paratively  limited  camera  for  output-recording,  one  with  an  ~2562  pixel  format.  The  repeated  pas¬ 
sage  of  the  data  through  the  system  is  required  by  the  need  to  move  the  small  camera  around  in  the 
output  plane  so  that  it  records  the  full  “image,”  i.e.,  the  Fourier  transform,  as  a  mosaic  of  smaller 
images.  (In  practice,  subsets  of  the  output  bandpass  are  electronically  selected  and  appear  within  the 
camera’s  stationary  format.)  This  compromise  has  been  accepted  because  the  goal  is  to  demonstrate 
the  principle  of  a  hybrid  architecture  in  the  most  expeditious  manner.  Although  performance  of  the 
prototype  will  be  suboptimal,  its  overall  speed  can,  in  principle,  be  easily  upgraded  by  using  a 
readout  camera  that  captures  the  entire  Fourier  spectrum  in  one  frame. 

The  architecture  described  below  is  realizable  in  a  variety  of  distinct  configurations,  each  of 
which  serves  a  different  purpose.  The  general  idea  of  a  2-D  frequency  format  has  existed  for  some 
time.  Mark  II  is  a  modification  of  a  design  that  was  initially  described  in  Refs.  10  and  1 1  but  never 
constructed.  The  reason  for  the  delay  between  conception  and  attempt  at  implementation  is,  essen¬ 
tially,  that  two  key  hardware  components  that  were  required  for  the  originally  intended  application  did 
not  exist.  These  components  were  an  output  camera  having  a  readout  speed  faster  than  any  that  exist 
even  today  (by  several  orders  of  magnitude)  and  a  chirp  device  capable  of  excursions  of  a  few  hun¬ 
dred  MHz/few  ms.  Thus  there  was  a  poor  match  between  available  technology  and  the  intended  pur¬ 
pose.  In  contrast,  for  the  intended  astronomical  applications,  Mark  II  is  a  high-performance  device 
that  is  available  with  current  technology.  Moreover,  the  fast  chirp  device  used  an  innovative  design. 
Having  incurred  the  primary  setup  cost  in  developing  Mark  I,  the  further  cost  for  Mark  II  has  been 
mininal.  We  now  turn  to  a  more  detailed  description  of  the  various  points  just  introduced. 

Expected  Performance 

The  expected  performance  of  Mark  II  derives  straightforwardly  from  the  parameters  of  key 
components.  Appendix  B  provides  the  predicted  figures  of  merit,  comparisons  with  the  Cray-2,  and 
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specifications.  The  camera  specification  (62500  pixels)  (Table  Bl)  shows  that  we  should  approach  a 
bandwidth  of  2 15  frequency  channels.  We  assume  that,  as  in  Mark  I,  the  limit  of  pixels/2  will  be 
attained.  Two  new  tellurium  dioxide  AOBDs  increase  the  center  operating  frequency  to  350  MHz 
(with  a  3  dB  bandwidth  of  200  MHz);  this  determines  the  maximum  signal  input  frequency  and, 
therefore,  the  data  rate.  The  dynamic  range  is  expected  to  be  -500,  which  is  a  substantial  factor  of 
20  increase  over  systems  that  are  of  a  fully  time- integrating  rather  than  hybrid  design.  This  results 
from  a  combination  of  the  low  light  bias  associated  with  the  space-integrating  dimension  and  the  use 
of  the  single-sideband  method  of  introducing  the  input  signal.  Thus  Mark  II  has  the  capability  of 
finding  weak  signals  in  the  presence  of  strong,  saturating  signals.  The  space-integrating  feature  also 
provides  a  gain  in  bandwidth  of  about  a  factor  of  four  over  purely  time-integrating  systems.  The 
throughput  (Table  B2),  which  results  from  considering  input  rate,  bit  accuracy,  and  the  number  of 
times  the  input  must  be  recycled  to  accommodate  the  output  camera,  is  about  double  that  of  a  Cray 
processing  the  same  transform. 

Design  and  Concept 

Like  Mark  I,  the  architecture  of  Mark  II  is  composed  of  two  separate  legs  of  an  interferometer; 
each  leg  produces  a  resultant  signal  at  the  output  camera  (Fig.  6).  In  a  wideband,  fast  processor  like 
Mark  n,  coherent  Fourier  integration  over  a  long  input  data  stream  is  desired.  The  wide  bandwidth 
is  accomplished  by  using  two  kinds  of  chirps.  In  Mark  II,  the  two  light  beams  carry  a  coarse  chirp 
(reference  leg)  and  the  input  data  (signal  leg),  both  introduced  by  using  the  350  MHz  AOBDs.  The 
coarse  chirp  leg  uses  a  fast  repeating  VCO,  which  is  the  most  demanding  chirp  generator  in  the  sys¬ 
tem.  The  device  must  sweep  100  MHz  ±  1  Hz  (we  plan  to  use  only  half  the  bandwidth  of  the  new 
AOBDs)  each  3  /is  and  return  to  the  same  starting  phase  ( ±  1  °)  each  repetition.  Without  return  to 
fiducial  phase  and  frequency  at  the  start  of  each  chirp,  the  integration  would  not  be  coherent.  The 
wideband  AOBDs,  together  with  the  spherical  lenses  that  follow  them  constitute  the  space-integrating 
portion  of  the  system.  The  narrow-band,  slow  chirps  are  introduced  into  both  legs,  propagating  in 
AOBDs  aligned  perpendicular  to  the  space-integrating  AOBDs.  The  geometry  and  operation  of  the 
horizontal  (slow  chirp)  dimension  is  similar  to  that  in  the  Mark  I  processor.  The  integration  is  car¬ 
ried  out  over  the  duration  of  the  slow  chirps.  The  interference  of  all  signals  produces  a  2-D  format 
Fourier  kernel  at  the  output  camera,  which  is  a  focal  plane  of  the  system.  The  effective  mathematical 
kernel,  however,  is  1-D:  the  processor  performs  a  continuous,  coherent  integration  over  only  one 
variable— time  (Fig.  7).  The  space-integrating  arm  produces  the  coarse-frequency  vertical  dependence 
in  the  output,  which  appears  as  a  succession  of  lines.  These  lines  are  associated  with  a  repeating 
delta  function  of  the  vertical  coordinate,  which  arises  from  recycling  the  fast  chirp.  The  horizontal 
direction  is  the  fine-frequency  direction  and  is  produced  by  the  slow  chirp  entered  into  both  arms.  A 
sinusoidal  input  will  produce  a  single  bright  line  (the  line  containing  the  input  frequency);  along  this 
line  will  be  a  single  bright  point  at  the  coordinate,  which  represents  a  fine  measure  of  the  input  signal 
frequency. 

Thus  four  AOBDs  are  in  the  design— two  accept  the  slow,  fine  chirp;  one  accepts  the  fast, 
coarse  chirp;  and  one  accepts  the  data  stream.  These  signals  are  folded  together  by  lenses  in  the  sys¬ 
tem;  the  result  can  be  represented  as  a  convolution  integral  [10,  11].  When  the  integral  is  evaluated 
to  calculate  expected  performance  at  the  output  plane,  the  result  is  an  expression  having  parts  that 
stand  in  one-to-one  correspondence  with  the  salient  features  of  the  output  format.  Appendix  C 
presents  a  complete  mathematical  representation  of  the  waveform  interactions  in  Mark  II  from  input  to 
output.  The  mathematical  treatment  was  verified  by  using  the  Symbolic  Manipulation  Program,  as 
described  of  the  end  of  this  section. 
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Fig.  6  —  Mark  D:  Two-dimensional  format,  hybrid  (space- 
and-time-integrating)  architecture.  Folded  one-dimensional 
Fourier  transform  processor,  a  signal  searcher  with  AO 
capability  for  recovery  of  smeared  signals. 
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Lens  performs  Fourier 
transform  (space  integration) 
of  fast,  repeating  chirp 
(5  ms  per  repetition, 

500  repetitions). 

FT  of  repeated  rect 
function  is  raster  of 
reference  frequency 
lines  with  range  of 
chirp  (300  to  400  MHz). 

Signal  beats  with 
one  frequency  on  one 
line,  forming  bright 
spot  on  bright  line. 


Time-integrating 
slow  chirps  (2.5  ms) 
are  imaged  to  the 
detector  plane, 
added  to  fast  chirp. 
Effective  range  is 
100  MHz/ 500  repetitions  = 
200  kHz.  Each  line 
receives  slow  chirp. 

Fast  chirp  is  reset 
to  same  phase  at 
beginning  of  each 
sweep;  time 
integration  is 
coherent  over  2.5  ms. 
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Fig.  7  —  How  Mark  II  performs  1-D  Fourier  transform  in  2-D 
(2-D  format  is  effectively  a  folded  1-D  Fourier  transform) 
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Coherence  Recovery  in  Two  Dimensions 

As  in  Mark  I,  implementing  the  CR  technique  in  Mark  II  requires  mixing  two  slow  unequal 
chirps  [6j.  The  2-D  format  does  not  affect  the  validity  of  the  method.  A  frequency-smeared  signal 
that  would  otherwise  have  overlapped  two  coarse  frequency  lines  is  recovered  at  one  point  on  one  line 
for  the  optimum  quadratic  transformation.  Moreover,  since  the  Mark  II  output  is  formatted  in  two 
independent  spatial  dimensions,  the  processor  can  be  designed  for  either  of  two  general  classes  of 
problems: 

•  applications  similar  to  the  astronomical  CR  problems  described  in  Section  1  where  the  2-D 
format  is  effectively  a  long  1-D  frequency  format  with  high  resolution  and  the  quadratic 
transformations  are  applied  consecutively;  or 

•  applications  (perhaps  not  astronomical)  in  which  rapid  CR  is  desired  and  low-frequency  reso¬ 
lution  will  suffice,  where  one  dimension  is  the  frequency  coordinate  and  the  other  simultane¬ 
ously  accommodates  parallel  processing  of  all  the  quadratic  transformations. 

The  next  two  sections  describe  architectures  of  these  types. 

High-Resolution,  Large  Bandwidth  Serial  Architecture 

In  the  high-resolution  Mark  II  processor,  the  CR  concept  can  be  realized  in  several  forms  that 
differ  in  how  the  slow  chirps  are  introduced.  The  slow  chirps  may  be  applied  in  at  least  two  addi¬ 
tional  geometrical  configurations  besides  the  “additive”  architecture  shown  in  Fig.  6  (see  Ref.  12  for 
a  discussion  of  additive  and  multiplicative  architectures).  First,  both  slow  chirps  may  be  introduced 
by  using  AOBDs  in  a  single  arm  of  the  architecture,  with  accompanying  lenses  rearranged  for  the 
required  imaging  and  spatial  Fourier  transformations.  The  effective  difference  in  acceleration,  which 
is  incurred  by  the  angular  compression  (from  frequency-dependent  Bragg  diffraction)  of  the  first  chirp 
when  traveling  the  distance  to  the  second  Bragg  cell,  can  be  compensated  by  a  base  difference  in 
acceleration  between  the  two  chirps.  Difference  in  acceleration  greater  than  this  frequency  offset  pro¬ 
vides  CR.  This  “multiplicative”  arrangement  is  apparently  not  described  in  the  literature.  An  alter¬ 
native  multiplicative  design  uses  one  slow  chirp  applied  through  an  AO  spot  modulator,  in  the  arm 
that  contains  the  slow  chirp  AOBD,  or  in  the  opposing  arm  [10,1 1].  In  any  one  of  these  three  config¬ 
urations,  the  difference  in  acceleration  between  the  two  slow  chirps,  which  constitutes  the  CR  portion 
of  the  chirp,  may  be  created  separately  as  a  third  chirp  and  single-sidebanded  with  the  signal  or  with 
the  coarse  chirp,  instead  of  appearing  in  one  of  the  slow  chirps.  The  latter  alternatives  for  CR 
emphasize  a  more  symmetrical  design  (with  both  slow  chirps  having  equivalent  acceleration),  but  they 
do  not  provide  a  materially  different  result.  Relative  advantages  of  each  arrangement  derive  primarily 
from  convenience  in  hardware  combinations  since  the  mathematical  descriptions  are  similar. 

In  all  of  the  foregoing  arrangements,  the  several  quadratic  transformations  required  to  achieve 
CR  are  performed  serially.  Thus,  for  103  transformations,  at  3  ms  per  run,  the  entire  set  of  transfor¬ 
mations  is  evaluated  in  -3  s.  The  merit  of  this  kind  of  processor  is  high  frequency  resolution  across 
a  large  bandwidth. 

Ultrafast  Parallel  CR  Architectures 

Mark  II  can  be  modified  to  perform  ~  103  time  transformations  in  parallel  with  a  reduced 
number  (—  103)  of  frequency  resolution  elements.  This  modification  would  provide  the  capability  of 
detecting  weak  signals  from  sources  undergoing  rapid  acceleration  or  from  signals  that  pass  through 
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an  inhomogeneous  medium  and  experience  an  attendant  frequency  smearing  that  is  roughly  linear  in 
time.  Since  all  ~  103  transformations  are  performed  simultaneously,  the  CR  integration  time  could  be 
as  fast  as  ~3  ms  with  the  current  design.  As  in  astronomical  applications,  such  a  processor  would 
require  a  speed  buffer  to  shift  the  original  low-frequency  signal  band  to  a  few  x  10  MHz,  100  MHz, 
or  GHz  bands,  depending  on  the  AOBD  bandwidth  of  choice.  Another  possible  application  for  the 
parallel  architecture  is  detection  of  weak  radio  signals  that  sweep  in  frequency  approximately  linearly. 
The  value  of  this  modified  version  comes  from  the  immense  speed  of  coherence  recovery  rather  than 
from  high  resolution  across  a  large  bandwidth  in  the  serial  design. 

The  modification  involves  electronically  offsetting  the  input  signal  by 
A/ signal  =  (i7AT)(A/fas, chirp  +  Max  {A/CR})  for  each  N  vertically  displaced  lines.  Hence,  the  signal 
is  simultaneously  diffracted  to  all  of  the  output  lines,  but  each  line  receives  a  different  quadratic 
transformation  (i.e.,  a  different  slow  chirp).  N  -  1  horizontal  lines  will  contain  the  signal  with  vari¬ 
ous  degrees  of  frequency  smearing  and  lower  intensities;  one  line  will  contain  the  signal  recovered  in 
a  single  channel.  The  fast  chirp  can  be  replaced  by  an  array  of  discrete  comparison  frequencies  that 
are  simultaneously  introduced  into  the  reference  leg;  this  choice  may  depend  on  ease  of  implementa¬ 
tion  within  the  bandwidth  of  interest.  To  minimize  intermodulation  products,  the  acoustic  modulation 
per  offset  must  be  kept  low.  The  same  geometrical  variations  in  slow  chirp  introduction  described  for 
the  high- resolution  design  are  possible  in  the  parallel  design. 

If  searches  for  (nearly)  linearly  frequency-swept  signals  are  performed,  the  family  of  chirp 
accelerations  can  then  be  adjusted  to  achieve  CR  on  frequency  drifts  of  perhaps  as  much  as  25%  to 
50%  of  the  signal  bandwidth.  The  essential  space- integrating  aspect  of  the  architecture  in  the  vertical 
dimension  is  retained,  while  the  integration  period  and  signal  band  are  easily  adjusted  for  specific 
applications. 

Symbolic  Manipulation  Program  as  Design  Tool 

Comparing  candidate  designs  is  a  basic  aspect  of  design  selection  in  developing  a  system  such  as 
Mark  13.  Examining  a  number  of  candidate  architectures  and  calculating  their  performance  provides  a 
variety  of  information.  The  calculation  overhead  is  considerable,  and  occasionally  the  architectures 
described  in  the  literature  contains  errors.  These  calculations  are  expedited  by  the  use  of  automated 
algebraic  processors. 

MACSYMA  and  SMP  are  two  examples  of  symbolic  manipulation  routines,  one  of  whose  main 
uses  is  the  checking  of  lengthy  and  tedious  calculations.  We  decided  to  develop  a  package  for  SMP 
that  would  calculate  and  confirm  expected  performance  of  AO  architectures.  Whereas  SMP  contains 
general  capabilities  for  evaluating  complex  expressions  by  using  rule-based  logic  and  identity 
libraries,  the  analysis  of  AO  architectures  called  for  developing  a  specialized  library  of  identities  and 
theorems  dealing  with  Fourier  transforms,  discontinuous  functions,  and  multiple  integrals. 

Beyond  performance  verification,  SMP  offers  other  advantages  worth  noting  and  exploring.  The 
SMP  calculation  constitutes  a  very  concise  algebraic  explanation  and  documentation  of  the  expected 
performance  of  the  system.  Certain  lines  of  the  calculation  correspond  to  the  operations  of  particular 
components.  Further  potential  also  exists  to  use  SMP  for  modeling  imperfect  performance.  The 
nominal  performance  of  a  component  may  be  replaced  by  a  representation  of  suboptimal  performance, 
and  the  result  may  be  propagated  through  to  later  stages.  SMP  facilitates  the  task  of  splitting  out  the 
state  of  the  optical  signals  generated  at  intermediate  steps  of  the  process.  Appendix  C  presents  an 
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SMP  calculation  that  simulates  Mark  II  performance,  symbolically  calculating  the  output  as  a  function 
of  the  input.  Results  of  selected  stages  are  shown  and  annotated  to  indicate  their  relationship  to  com¬ 
ponents  of  the  processor. 

FUTURE  DEVELOPMENTS 

The  program  orientation  remains  development  of  processors  for  astronomical  data,  either  on  the 
ground  or  in  orbit.  At  some  future  point  it  may  become  important  to  explore  packaging  of  such  pro¬ 
cessors  for  space  flight. 

An  ambitious  extension  of  Mark  II  is  contemplated  that  could  provide  precision  equal  to  that  of 
a  digital  machine.  We  are  considering  developing  a  digital  acousto-optical  array  processor  in  which 
some  aspects  of  the  design  resemble  a  generalized  Mark  II  architecture.  This  array  processor  would 
be  capable  of  linear  operations  with  throughputs  (approaching  1012  bits/s)  almost  two  orders  of 
magnitude  greater  than  those  of  state-of-the-art  digital  supercomputers;  it  would  overcome  the  limita¬ 
tion  of  inadequate  precision  that  is  inherent  in  many  analog  processor  systems. 
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Appendix  A 

COHERENCE  RECOVERY  OF  A  FREQUENCY-MODULATED  SIGNAL  BY 
USING  QUADRATIC  TIME  TRANSFORMATIONS 


For  an  emitting  source  undergoing  arbitrary  acceleration,  a  grid  of  transformations  must  be 
applied  in  a  space  of  arbitrary  dimension  to  achieve  optimal  CR  in  only  one  transformation  of  the  set. 
For  many  applications,  a  quadratic  approximation  over  a  short  portion  of  the  trajectory  of  points  in 
/i-parameter  space  may  be  sufficient  to  detect  an  otherwise  unobservable  signal.  The  great  value  of 
this  method  is  that  the  quadratic  approximation  requires  only  one  parameter,  thereby  resulting  in  a 
substantial  saving  of  computation  time  and  resources.  In  most  cases  of  astrophysical  interest,  the  sig¬ 
nal  arises  from  a  source  in  a  binary  orbit,  which  is  an  ellipse.  A  circular  approximation  to  an  ellipti¬ 
cal  orbit  provides  an  adequate  example  of  the  contraction  of  parameter  space  from  three  variables  to 
one  and  does  not  detract  from  the  validity  of  the  treatment. 

Coherence  Recovery  in  Binary  Stellar  Orbits 

For  a  circular  orbit  the  observed  signal  is  advanced  or  retarded  by  a  time  of  maximum  ampli¬ 
tude  A,  the  projected  radius  of  the  orbit  expressed  in  light  seconds,  with  respect  to  the  center-of-mass 
frame.  The  system  revolves  with  frequency  Q  =  2x7 Porb,  where  Porb  is  the  orbital  period.  The 
phase  lead/lag  may  be  expanded  in  a  Taylor  series  to  second  order  (an  arbitrary  orbital  phase  4>o  is 
neglected  for  simplicity): 

6<t>  =  A  •  sin  (Qt)/Parb  =  ,4  [sin  (0  T0)  +  Qcos  (QTa)t  -  Q2  sin  (SlT0)t2 /2]/Polh.  (Al) 

T„,  the  time  about  which  the  function  is  expanded,  can  be  taken  to  be  the  center  of  observation.  The 
first  term  of  the  expansion  represents  this  fixed  point  in  phase  with  respect  to  the  orbit;  it  is  of 
interest  only  in  absolute  phase  studies.  The  second  term  represents  an  average  Doppler  shift  for  the 
interval  of  observation  and  contains  no  broadening  effect  since  it  is  linear  in  time;  pulses  arriving 
P puise  seconds  apart  will  be  compressed  or  expanded  by  ±  [1  +  Q/i  •  cos  (OTo)],  which  is  a  constant 
factor.  The  last  term  constitutes  the  lowest  order  smearing  or  broadening  effect  in  phase  of  arrival 
and,  therefore,  in  frequency.  This  term  is  quadratic  in  phase  and  linear  in  frequency,  and  it  attains  a 
maximum  when  the  argument  of  sin  is  r/2.  Hence  when  the  source  exhibits  zero  Doppler  shift 
(second  term),  the  broadening  is  greatest  and  can  be  of  positive  or  negative  “curvature.”  A  qua¬ 
dratic  approximation  applied  to  the  time  of  arrival  of  the  sampled  events  will  compensate  for  the 
phase  lead  or  lag: 


t'  =  t  +  0 2 A  •  sin(flro)t2/2  =  t  +  at2 . 


(A2) 
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The  number  of  transformations  required  to  cover  positive  or  negative  curvatures  for  observations  of 
worst  case  (Q T0  ~  ir/2)  is 


^  tran  =  <W/6a  =  (AQ2/2)/(PCTil/2T2).  (A3) 

The  expression  for  6a,  derived  below,  assumes  a  program  that  searches  for  optimal  recovery,  which 
is  desirable  when  the  signal  may  be  barely  detectable.  Equation  (A3)  can  be  rewritten  as 

Ntm  =  A(2irT/P  orb)2  IP  pu)se  (A4) 

if  P puise  *s  near  Pcnt  =  P Nyq.  Otherwise,  the  required  number  of  transformations  is  reduced  by  a 
factor  P puise/^Nyq.  In  a  search  problem,  Ppulse  is  not  known  a  priori,  and  therefore  the  maximum 
number  of  transformations  is  required.  Note  that  Nusn  is  proportional  to  (2-kT/P orh)2 .  Moreover, 
2irT/Polh  must  be  small  for  the  approximation  to  be  valid.  As  2itT/Porb  grows,  the  higher  order 
terms  become  more  significant,  and  quadratic  CR  becomes  less  optimal.  The  third-order  term  is 
important  when,  assuming  a  phase  tlT0  =  n/4, 

A&T2/ 6  -  AQ2T2/2, 

therefore  requiring  the  constraint  that  T/P orb  «  3/2 ir  for  the  quadratic  approximation  to  still  be 
useful.  Note  that  when  the  quadratic  effect  approaches  minimum,  the  third  order  begins  to  dominate. 
Conversely,  near  the  region  of  maximum  quadratic  effect,  the  third  order  is  near  minimum.  Hence 
the  most  interference  with  quadratic  recovery  from  the  third-order  term  does  arise  near  phase  ir/4. 
The  restriction  on  practical  integration  times  can  be  recast  into  a  constraint  on  orbital  periods 
searched.  Fortunately,  a  large  portion  of  parameter  space  is  made  accessible  by  quadratic  transforma¬ 
tions,  extending  useful  integration  times  in  binary  pulsar  searches  by  one  to  three  orders  of  magni¬ 
tude. 


We  now  calculate  the  critical  5a  required  to  recover  a  broadening  of  the  signal  by  precisely  one 
channel  for  a  period  Pcrit.  Taking  the  derivative  of  the  quadratic  transformations  (Eq.  A2),  we  find 

dt'  =  dt(  1  +  26af)  -  6a  =  (dt'/dt  -  l)2f. 

Evaluate  for  the  total  integration  time  T  and  require  the  transformation  to  stretch*  the  time  series  by 
the  amount  Pcnt: 

6a  =  [(Pcrit  +  T)/T  -  11/2 T  =  Pcri(/2r2  .  (A5) 

♦Incrementing  the  integration  time  by 

4T  =  P Nyq  ■  1  cycle  at  the  Nyquist  frequency 

adds  one  channel  in  the  resulting  Fourier  power  spectrum.  Let  Npt  =  77 2 r,  =  number  of  channels  in  power  spectrum  for 
integration  time  T,  where  r,  is  the  sample  interval.  Then 

T‘  =  T  +  2r,  =  T  +  P Nyq,  and  therefore  Np, '  =  Np,  +  1 . 

Transformation  (T  —  T1)  of  a  time  series  from  one  that  contains  a  coherent  signal  to  one  where  the  signal  begins  at  PNyst ' 
and  ends  at  P Ny9'  results  in  broadening  by  exactly  one  channel.  The  distribution  of  spectral  power  within  the  two  channels 
depends  on  the  specific  transformation. 
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Coherence  Recovery  Realization  In  Time-Integrating  Acousto-Optical  Fourier  Processors 

In  a  time-integrating  AO  Fourier  processor,  two  chirps  of  equivalent  acceleration  and  duration 
modulate  a  coherent  light  source  and  are  heterodyned  to  create  a  Fourier  kernel.  An  additional  phase 
with  quadratic  time  dependence  is  manifested  when  the  chirps  have  slightly  different  acceleration. 
The  modified  kernel  can  effect  exact  CR  of  a  signal  in  a  time  series  that  has  been  quadratically  time 
transformed;  approximate  recovery  is  realized  for  other  functions.  For  signal  5(r),  chirps  Cx  and 
Ct*  (complex  conjugates),  delay  t,  chirp  start  frequency  /0,  and  chirp  accelerations  a  =  a', 

S(v)  =  f  S(t)Cx(t,T)Cx*(t,  r)dt 

=  j  S(t)  exp  j[(t  +  T){a(t  +  t)  +  2ir/o)]  exp  j[-(t  -  T)\a'{t  -  r)  +  2 irf0}]  dt 

=  j  S(t)  exp y[47r/0T  +  2(a  +  a^tr  +  (a  -  a ^(t2  +  r2)]  dt.  (C6) 

Letting  a'  =  a  —  5a  and  a"  =  ( a '  +  a)/2  yields 

S(v)  =  (  S(t)  exp 7[4tt/ot  4-  4 a"tr  +  Sa(t2  +  r2)]  dt.  (C7) 

J  (1)  (2)  (3) 

The  unequal  chirps  result  in  the  terms  labeled  in  Eq.  (C7):  (1)  the  usual  acoustic  fringe  carrier 

(time-independent);  (2)  a  Fourier  kernel  slightly  modified  from  the  equal  chirp  case,  4a tr  —  4a "tr; 

and  (3)  an  additional  quadratic  term  in  t.  The  r2  dependence  in  the  third  term  represents  an 
irrelevant,  fixed  spatial  fringe  pattern  that,  in  practical  realizations,  is  of  vanishingly  small  magnitude 
compared  to  terms  of  interest.  The  t2  phase  in  the  integrand  can  remove  a  -t2  phase  in  a  signal. 
Let  the  signal 

5(f)  =  exp  (J2irvt)  —  Sit)  =  exp  j{2vvt  +  yt2).  (C8) 

When  5a  =  -y,  coherence  is  recovered  exacdy,  and  Eq.  (C7)  becomes 

5(v)  =  exp  j[4irf0T]  J  5(0  exp  ;[4a"fr]  dt,  (C9) 

which  is  the  usual  time-integrating  Fourier  AO  formulation,  where  the  relationship  between  AOBD 
delay,  r,  and  frequency  is  v  =  2a"T/r.  In  many  cases  the  quadratic  transformation  will  be  the  best 
second-order  approximation,  and  certainly  the  most  straightforwardly  generated  analog  implementa¬ 
tion  of  CR,  over  a  short  interval  to  a  more  complex  function,  e.g.,  FM  modulation  arising  from 
motion  over  a  portion  of  an  orbit  or  trajectory. 
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Appendix  B 

MARK  H:  FIGURES  OF  MERIT  AND  SPECIFICATIONS 


Table  B1  —  Mark  II  Figures  of  Merit 


Frequency  resolution  elements 

—  215  (62500  pixels  =  prese-’  camera)/2 

Input  data  rate 

-350  MHz 

Dynamic  range 

-500 

Signal  to  noise 

-  25  (higher  with  increased  light  level) 

Dynamic  range  increase  (factor  >  20)  over  two-dimensional  time-integrating  systems  results  from  use 
of  space  integration  in  coarse  frequency  dimension;  this  provides  the  capability  to  detect  very  weak 
signals  in  the  presence  of  strong,  saturating  signals. 

Data  rate  increase  over  Mark  I  results  from  use  of  Te02  Bragg  cells  with  center  frequencies  of  350 
MHz  to  introduce  signal  and  coarse  reference  chirp. 

Space  integration  also  enables  gain  of  -4  x  over  two-dimensional,  purely  time-integrating  systems 
in  useful  input  bandwidth  of  Te02  cell. 

System  throughput  (bits/s),  apart  from  postprocessing,  is  approximately  two  equivalent  Cray-2s,  by 
throughput  comparison  on  FFT  algorithm.  The  strictly  acousto-optic  portion  of  the  system  performs 
the  transform  —500  times  faster  than  Cray-2,  but  with  less  precision  (see  Table  B2). 


Table  B2  —  Throughput  Comparison  with 
State-of-the-Art  Electronic  Computer 


Cray-2: 

64-bit  accuracy 

106 -point  FFT  in  1.24  s 

—  (106 /2-point  spectrum) 

—  27  mbits/s  throughput 

AO  Processor: 

estimate  -  4-bit  accuracy  (initially) 

-  106  points  in  -2.5  ms  with  present  camera 
and  2  pixels  per  frequency  channel 

—  1/16  of  spectrum  recorded 

-  1  Gbits  (=4  bits  x  0.35  GHz)  input/16 
=  67  Mbits/s  throughput 

-  2  equivalent  Crays  (working  speed) 

♦Accuracy  increases  with  investment  in  improving  light  budget,  postprocessing. 
tAO  processor  does  transform  -500  times  faster  than  Cray-2. 
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Table  B3  —  Mark  II  Special  Components 

2  Specially  engineered  Te02  Acousto-optic  Deflectors 
200  MHz  @  3  dB,  5#ts  aperture 
Model  4350-1  Crystal  Technology 

1  Input  Signal  Driver  for  4350-1  Deflector 
Model  1350-1  Crystal  Technology 

1  Fast  chirp  VCO  2.5-/ts  sweep,  0.5-jts  dead-time 

phase-locked  to  1  ° ,  <  ±1  Hz  @  each  sweep  start 
100  MHz  sweep  uses  center  half  of  4350-1  BW 
SO-350-200-XX  Andersen  Laboratories 

(Unit  will  be  —  1.5  to  3  orders  of  magnitude 
faster  than  available  digital  sweeps} 
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Appendix  C 

VERIFICATION  OF  THE  MARK  H  DESIGN  BY 
USING  THE  SYMBOLIC  MANIPULATION  PROGRAM 

We  performed  an  automated  verification  of  the  Mark  II  acousto-optic  design  by  using  the  Sym¬ 
bolic  Manipulation  Program  (SMP).  These  computations  served  two  purposes.  First,  it  was  neces¬ 
sary  to  verify  that  the  processor  would  function  according  to  specifications  after  modifying  the  origi¬ 
nal  Bader  [Cl]  design.  Minor  errors  or  misprints  in  Ref.  Cl  were  uncovered  in  this  process.  These 
errors  result  in  minor  adjustments  (phase  shifts)  to  the  mathematical  representation  and  turn  out  to  be 
irrelevant  to  the  functioning  of  the  processor.  The  corrected  versions  were  verified  by  the  SMP  com¬ 
putations.  Second,  we  were  interested  in  exploring  further  modifications  to  the  design,  e.g.,  the  reali¬ 
zation  of  coherence  recovery  (CR)  in  a  two-dimensional  processor.  SMP  provides  the  capability  of 
perturbing  functions  in  a  previously  defined  stream  of  calculations  at  any  stage,  thereby  obviating 
manual  reanalysis  of  the  mathematical  problem  and  simplifying  evaluation  of  design  changes. 

An  interpretation  of  the  SMP  input  and  output  for  the  Mark  II  design  follows.  Equations  that 
have  one-to-one  correspondence  with  “operations”  are  numbered.  These  are  not  necessarily  in  one- 
to-one  correspondence  with  hardware  elements  since  such  elements  may  introduce  more  than  one  fun¬ 
damental  transformation.  Other  non-numbered  equations  are  given  for  clarity. 

The  numbering  system  for  the  following  equations  divides  them  into  four  groups  according  to 
how  they  function  in  the  calculation.  Equations  (5,)  are  associated  with  inputs  in  the  “signal”  leg  of 
the  processor;  Eqs.  (C,)  characterize  chirps  introduced  in  the  other  (“reference”)  leg  to  analyze  the 
input.  Mathematical  identities  used  to  compute  the  output  are  not  generally  numbered.  Finally,  Eqs. 
(/,)  describe  the  results  produced  by  the  processor,  corresponding  to  convolution  integrals  in  the 
detector  or  output  leg  of  the  processor.  (See  Fig.  6  for  the  system  overview  showing  the  three  legs.) 

In  one  leg  of  an  interferometer,  the  signal  S(t)  is  introduced  by  using  an  acousto-optic  beam 
deflector  (AOBD)  of  aperture  D.  S(t)  represents  the  amplitude  of  the  acoustic  information  (i.e.,  the 
data— for  simplicity,  one  unmodulated  rf  carrier  within  the  Bragg  cell  [AOBD]  bandwidth).  The  trav¬ 
eling  acoustic  wave  is  represented  by  S(t  -  £/v),  where  £  is  the  deflector  coordinate  (within  the 
cell),  approximately  transverse  to  the  light  propagation  vector,  and  v  is  the  acoustic  velocity.  The  dif¬ 
fracted  light  amplitude  A(x ,  t)  is  proportional  to  the  acoustic  signal  amplitude. 

Signal  Input  Leg: 

The  amplitude  spectrum  U(f)  of  S(t)  is,  by  definition, 

(/(/)  =  j  dt  S(t)  e\p\i2irft\. 

(acoustic  signal  propagates  through  AOBD) 

A(S,t)  =  S(t  -  i/v)  (SI) 
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(By  Bragg  diffraction,  acoustic  signal  impressed  onto  light  beam  which  passes  through  AOBD  of 
aperture  D;  signal  is  now  light  amplitude) 

A(£,  t)  =  S(t  -  f/v)  rect  (f/D), 

where 

rect  (f/D)  =1  for  -D/2  <  f  <  D/2 

=  0  otherwise.  (S2) 

(light  signal  Fourier  transformed  by  lens) 

A(x,  t)  =  j  df  S(t  -  f/v)  rect  (f/D)  exp{ -2xax£},  (S3) 

where  x  is  the  vertical  (course  frequency)  detector  coordinate,  a  =  1  /(KF),  X  =  light  wavelength, 

and  F  =  focal  length  of  the  Fourier  transform  lens.  Equation  S3  assumes  that  imaging  from  the  sig¬ 
nal  AOBD’s  Fourier  transform  plane  to  the  detector  plane  is  one-to-one.  A(x,  r)  is  in  the  detector 
plane,  which  is  the  imaged  Fourier  transform  plane  of  the  signal  AOBD. 

In  the  following,  recognize  the  auxiliary  variable  z  =  f/v  ( z  is  a  scaled  Fourier  transform  coor¬ 
dinate  of  the  time  t).  Equation  (S3)  is  reduced  to  another  form,  and  an  example  of  the  convolution 

theorem  is  explicitly  derived.  First,  let  t  —  r  -  £/v  (using  shift  theorem), 

U(f)  exp{i'2rr£/v)  =  \  dt  S{t  -  £/v)  exp(/2ir/r), 

and  taking  the  Fourier  transform  of  both  sides, 

f  df  U(f)  exp(-/27r/(r  -  £/v)|  =  f  df  U(f)  exp(-/'2ir/r|  exp{/2ir/£/v|  =  S(t  -  f/v). 


Substituting  this  result  into  Eq.  (S3)  gives 

(Shift  theorem  used  to  transform  time  shift  to  phase  shift) 

A(x,  t)  =  j  d(  |  df  U(f)  exp[-i2irft]  exp{/'2Tr/f/v)  rect  (f/D)  exp| -/2iraxfl .  (S4) 

Now  remembering  that  the  FT  of  rect(jc)  is  sincfxAO  =  sin  (rX)/irX,  df  =  v  •  dz,  and  recognizing 
that  a  factor  1/v  arises  in  the  transformation  (/(/)—  U(zv)/v, 

A(x ,  0  =  J  v  •  dz  U{zv)/v  exp| -i2irzvr|  j  df  rect  (f/D)  exp( -/2irf(ax  -  z)), 

(Integrate  out  f  dependence) 

A( x,  t)  =  j  dz  exp{ -Hirzvt |  U(zv)  sinc[irD(a£Jt  -  z)].  (S5) 
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Recapitulating,  A(x,  t)  is  the  diffracted  light  amplitude  containing  the  information  in  the  signal 
at  the  detector  plane,  having  passed  through  the  AOBD  of  width  D  and  having  been  Fourier 
transformed  by  a  lens.  One  effect  of  the  latter  is  represented  by  the  phase  factor  ax  in  the  sine  func¬ 
tion;  also  U(zv)  is  the  optical  Fourier  transform  of  the  input  signal. 

Reference  Leg; 

In  the  other  leg  of  the  interferometer,  a  second  AOBD,  having  properties  identical  to  the  signal 
AOBD,  is  driven  by  a  periodic  chirp  of  period  T , .  The  chirp  has  bandwidth  A /,  and  center  fre¬ 
quency  fc,  parameters  that  span  the  potential  range  of  signal  frequencies.  This  coarse  chirp  is 
represented  by 

(Acoustic  chirp  signal  produced  by  chirp  generator) 

C\(t)  =  £  rect  [(/  -  nTx)/T\  \  exp{i[a(f  -  nTxf  +  b(t  -  nr,)]),  (Cl) 


where  a  =  irA f\/T\  and  b  =  2ir/min,  and  a  constant  phase  term  <t>n  in  the  b  coefficient  is  suppressed 
for  convenience  (notice  that  the  coarse  chirp  generator  must  supply  chirps  each  beginning  with  this 
phase,  <t>„  =  <Aconst ) .  This  repeating  chirp  propagates  through  the  AOBD  and  is  Fourier  transformed 
by  a  lens  identical  to  the  lens  that  transforms  the  signal  beam.  We  proceed  to  determine  the  ampli¬ 
tude  of  the  chirp  light  field  in  the  detector  plane: 

(Chirp  propagates  through  AOBD  as  acoustic  signal) 

C,(r,  £)  =  £  rect  [(t  -  nT x  -  £/v)/r,] 
n 


■  exp{i[a(f  -  nT\  -  $/v)2  +  b(tnT ,  -  £/v)]).  (C2) 

(By  Bragg  diffraction,  chirp  impressed  onto  light  beam  that  passes  through  AOBD  of  aperture  D; 
chirp  is  now  carried  by  light  amplitude) 

r(t,  £)  =  £  rect  [(t  -  nTx  -  ^/v)/?,]  •  rect  (£/D) 

n 

■  expj/[a(f  -  nT{  -  $/v)2  +  b(t  -  nT}  -  £v)]|.  (C3) 

(The  optical  Fourier  transform  of  chirp  is  performed  and  imaged  one-to-one  to  the  detector  plane) 

R(t,  x)  =  J  d$  r(t,  0  exp{  -/ 2ira^x) 

=  J  d£  exp(-/2ira£x)  £  rect  [(/  -  nT\  -  ^/v)/7' , ]  •  rect($/D)  (C4) 

n 

■  exp{i[a(f  -  nT\  -  £/v)2  +  b(t  -  nT \  -  £/v)]). 
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In  correspondence  with  the  Convolution  Theorem,  make  the  associations 

f  =  Y  rectt(*  -  nTi  -  Z/v)/Tx]  exp{i[a(r  -  nTx  -  £/v)2  +  b(l  -  nTx  -  £/v)]} 
and  g  =  rect(£/D),  and  the  following  result  is  straightforwardly  obtained: 

(Use  Convolution  Theorem  to  rearrange  chirp) 

R(t,  jc)  =  D  |  dx '  V(x  -  x ',  t)  sinc[7ra£)jc  'J  , 

where 

V(x,  t )  =  di,  exp{/[a(r  -  nTx  -  £/v)2  +  b(t  -  nTx  -  ^/v)  -  2ira£x]} 

•  rect[(r  —  nTx  —  S/v)/T1],  (C5) 

In  K(jc,  t)  make  the  substitution  z  =  t  -  nTx  -  £/v,  -vdz  =  d£  (note:  rect  function  is  symmetric 
and  therefore  introduces  a  minus  sign,  absorbing  the  minus  that  arises  from  —vdz): 

(Variable  change  simplifies  notation) 

V(x,  t)  =  v  Y,  exp{-2?ra vx(t  -  nTx)} 

n 

•  j  dz  expjiaz2  +  iz[2iravx  +  b\\  rect[z/T,].  (C6) 

Again  using  the  Convolution  Theorem,  transform  the  integral  in  Eq.  (C6).  The  problem  can  be 
formulated  such  that  the  dual  variables  of  the  Fourier  integral  are  z  and  avx  +  b  Hr.  In  correspon¬ 
dence  with  the  convolution  theorem,  make  the  associations 

/  =  rectlz/T!)  and  g  =  exp{/oz2). 

The  phase  factor  izb  will  appear  in  the  transform  of  /  as  a  shift  and  in  the  transform  of  g  as  a  phase 
term. 

To  see  how  to  arrive  at  Eq.  (C7),  first  obtain  the  Fourier  transform  of  g  alone,  then  generate  the 
convolution  of  the  transforms  of  /  and  g.  The  trick  is  to  use  the  transform  property  of  a  Gaussian. 

Let  avx  +  b  Hr  =  u.  Then, 

J  dz  expfiaz2}  exp[i2ruz)  =  j  dz  expf-x  [~iaz2/r  -  2 iuz]). 
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Completing  the  square, 

j  dz  exp{-x  [-iaz2 /r  -  2 iuz  +  ( r/ia)u 2  -  (r/ia)u2]}  = 
exp(-/(7rw)2 /a\  |  dz  exp{-7r[z(  — ia/r)l/2  -  iu/(-ia/r)v2]2}  = 

Letting  z'  —  z(-ia/r)l/2  -  iu/(-ia/r)l/1 ,  dz  =  (— ia/r)~l/2dz we  get 

exp{-/(7TM)2 /a)  j  (iTr/a)l/2dz'  expf-xz'2}  =  ( ir/a)l/ 2  exp{-/(T«)2/a} . 

Now  transform  /  and  multiply  with  the  folded  transform  of  g,  include  the  shifts  arising  from  exp (izb, 
substitute  avx  +  b  /2r  for  u,  and  arrive  at  a  new  expression  for  the  integral  in  Eq.  (C6): 

(Identity  to  facilitate  use  of  Convolution  Theorem  in  Eq.  (C6)) 

T,  K  j  dx'  exp{-i(xav)2(x  -  x'  +  b  /2rav)2 /a]  •  sincIx-T]  (vox  ’  +  b/2r)\,  (C7) 

where  K  =  ( ir/a)xn .  Now  proceed  to  simplify  Eq.  (C7)  by  an  approximation,  following  Bader 
exactly.  Expanding  the  exponent  in  Eq.  (C7), 

exp{-/(27rav)2[x2  +  x'2  +  b2/(2rav)2  +  xb/vav  -  2xx'  -  x  'b/rav]/Aa] 

The  (x r)2  term  is  negligible  in  the  region  over  which  the  sine  function  is  significant:  at  one  Rayleigh 
distance  x '  =  l/(vT{a),  and  therefore 

(r avx')2 /a  =  ir/fA/iTj)  =  r/Nx  «  1, 

where  we  have  used  the  definition  a  =  rAf\/T{  and  where  is  the  number  of  resolvable  frequen¬ 
cies  in  the  coarse  chirp.  Moving  the  non-*'  terms  outside  the  u  =  x  4-  b/2rav,  Eq.  C7  looks  like 

T\  K  expl-i*}  j  dx’  tx\>[i2rx’u\{2vav)2 /Ara})  sincTTj  [vox '  +  b/2r]. 

Now  letting 

x  "  =  vox '  +  b/2r,  x '  =  x"/va  -  b/2rva,  dx'  =  dx  "/va,  we  get 
T\  K/va  exp(-i%(  j  dx  "  exp{i2iru(x "/va  -  b /2rva)(2rva)2 /Ara]  sinclirTjx"]. 

Now  let 

w  =  x  "(2rav)2  /Arava  =  x"rva/a,  dx"  =  ( a/rva)dw : 

T |  Ka/r(v a)2  exp{-/^|  exp{ -irbavu/a]  J  dw  exp(j2irwu)  sinc[xriaw/Tva]. 
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Performing  the  Fourier  integral  and  inserting  u  =  x  +  b/2nva  results  in 

K/nva  exp{ -ib2/2a]  exp{  —  z cxp[—inbavx /a]  rect  [Tav(x  +  b/2nav)/Txa] . 

Let  c  -  (in /a)1'2  ■  [a /net]  exp  {-  ib2 /2a].  Then  we  have  from  above  substitutions  and  Eq.  (C6), 

V(x,  t)  =  c  exp( -i[27ravx  +  h]2/4a)  exp[-inbavx/a]  rect  [(7ravx  +  b/2)/Txa] 


•  £  exp[  —i 2navx(t  —  nTx)]. 

n 


Now,  using  AfxTx  =  Nx  =  ctD2 ,  Tx  =  D/v,  a  =  nA fx/Tx,  and  b  =  2nfmin,  define  the  rect 
argument  with  respect  to  /min  (set  r  =  f^/  Af\ ;  note  that  Bader  defines  r  =  fcen/Afx;  either  way, 
the  limits  are  consistent,  see  Eq.  (Cll)): 

rect  [(7ravx  +  b/2]/Txa]  —  rect  \{nav2x/a  +  bv/2a]/D]. 

Since  7rav2/n  =  n(AfxTx/D2)  (D2 /f\)/(nAfx/Tx)  =  1,  bv /2a  =  Dfmm/Afmm/Afx  =  rD,  then 
we  have 

(Several  transformations  to  evaluate  integral  from  Eq.  (C6)) 

V(x,  t)  =  c  exp{-/[2iravx  +  b]2/4a]  exp{-i7r&ciivx/a)  rect  [(x  +  rD)/D ]  (C8) 

•  £  exp{-i2navx(t  -  nTi)}, 


where  an  extra  phase  (exp[inbavx  /a])  is  discovered  that  does  not  appear  in  Bader’s  derivation.  The 
phase  is  physically  meaningful:  it  represents  the  discontinuity  from  one  coarse  frequency  line  to  the 
next  and  arises  from  a  combination  of  the  lens  Fourier  transform  and  the  fast  chirp.  It  is  inconse¬ 
quential,  since  we  desire  only  the  power  and  not  the  amplitude  spectrum. 

Continuing,  use  the  identity  £  exp(ian)  =  2n  £  6(a  —  2np) , 

n  p 

where  6(x)  is  the  Dirac  delta  function.  Make  the  association  2navxTx  =  a  to  arrive  at 
(Convert  sum  of  exponentials  to  [approximate]  sum  of  5s) 

K(x,  t)  =  c  expj -i[2iravx  +  b]2/4a]  exp{-inbavx/a]  rect  [(x  +  rD)/D] 


•  exp[-i'2Tavxr)  £  5(x  -  p/avTx). 

p 


(C9) 
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Now  combine  Eqs.  (C 9)  and  (C5),  setting  C  =  cD: 

R(t,  x )  =  C  j  dx '  exp{-i[2xav(x  —  jc1)  +  b]2  /4a}exp[-iirbav(x  -  x^/a] 


•  reel  {[(*  —  x')  +  rD]/D)  expj -i2votv(x  -  x')/}  sinc[xaD;r1 
■  Y  6[(r  -  x')  -  p/avTx]. 

The  delta  function  replaces  all  (x  —  x')  with  p/avT i  and  the  x 1  in  the  sine  argument  with 
x  -  p  /avT  |  to  yield 

R(t,  x)  =  C  Y,  exp [-ilZ-irp/Ty  +  b]2/4a)  exp[-iirbp /T{/a\  exp(-i2Trpr/7',} 

p 

•  rect  {[(/avT1!  +  rD]/D }  sine  [iraD{x  -  p/avT\)]. 


Rearranging  and  noting  again  that  vTx  =  D  and  N \  =  aD2,  we  obtain 
(Use  delta  function  and  identities) 

R(t,x )  =  C  Y  exp[-i2vpt/Tx  -  -  i[2vp/T\  +  b]2/4a\  txp[-irbp /aTx j 

p 

•  xecl\p/Nx  +  r]  sinc[7r7'1(avx  -  p/Tx )].  (CIO) 


The  effect  of  the  (stationary)  quadratic  phase  term,  exp{-/[2xp/7'i  +  b]2/4aj,  is  negligible 
and  irrelevant  if  a  power  spectrum  is  desired.  Also,  the  extra  phase  term  discovered  in  Eq.  (C8)  is 
similarly  irrelevant.  If  an  amplitude  spectrum  with  correct  phases  is  required,  the  effects  of  these 
terms  can  be  compensated  by  a  specially  designed  lens.  The  stationary  phase  terms  are  suppressed 
until  the  end  of  the  derivation.  The  simplified  form  results: 

(Neglecting  stationary  phases) 

R(t,  x )  =  C  Y  exp(-2irpt/Tx )  rect  \p/N ,  +  r]  sinclxT] (avx  -  p/Tx)),  (Cl  1) 

p 

where 


Nxr  <  p  <  Nx(r  +  i),  and  r  =  /min/A/,; 


equivalently. 


/min  —P /T X  ~  /min  A/l  • 

This  is  Bader’s  Eq.  (7),  the  curvature-corrected,  coarse  chirp  reference  wave  at  the  detector  plane. 


29 


WOOD,  NORRIS,  AND  SMATHERS,  JR. 


In  our  modification,  the  slow  chirp  multiplies  both  legs  (with  opposite  directions  of  propagation) 
by  using  AOBDs.  The  resulting  deflections  are  in  the  y  coordinate.  The  slow  chirps  propagating  in 
the  AOBDs  with  acceleration  /3  (minus  sign  for  signal  leg,  plus  sign  for  reference  leg)  and  start  fre¬ 
quency  /0  are 


C2  =  exp {/(/  ±  T)(P(t  ±  t)  +  27r/o)}. 


Detector  Output  Leg: 

Applying  the  slow  chirp  to  the  signal  leg  (Eq.  (S5))  yields  the  time-dependent  signal  contribution 
to  the  light  amplitude  at  the  detector: 

(Signal  amplitude  multiplied  by  slow  chirp) 

A(x,  t)  -  exp{i(r  -  r)[/3(?  -  r)  +  27r/0]}  j  dz  expf-j2irzvr}  U(zv)  sinc[7r£>(ax  -  2)].  (II) 


Similarly,  the  slow  chirp  multiplying  the  coarse  chirp  reference  (Eq.  Cll)  yields 
(Coarse  chirp  reference  multiplied  by  slow  chirp) 

R(t,  x)  =  C  £  cxp{  - ilirpt /Ti  \  rect  \p/N\  +  r]  sinclftTi (crar  —  p/T\)]  (12) 

p 

■  exp{  i(t  +  r)(0(t  +  t)  +  2ir/0)j. 

The  portion  of  the  light  bias  offset  arising  from  the  signal  leg  of  the  interferometer  is  determined  by 
taking  the  real  part  of  the  square  of  the  light  amplitude  from  Eq.  (II),  integrating  for  the  duration  of 
the  slow  chirp  T2,  and  making  the  substitution  z  ~  f/v: 

(Signal  leg  contribution  to  the  light  bias) 

Ei  =  T2  J  df  sinc[7rD(ax  -  f/v)]2  U(f).  (13) 


A  similar  calculation  for  the  light  bias  arising  from  the  coarse  chirp  leg  yields 
(Coarse  chirp  leg  contribution  to  the  light  bias) 

E2  =  C  T2  sinc[ir(avx7'1  -  p)]2  .  (14) 

Within  constants  (which  are  adjustable  by  varying  the  light  levels  in  the  two  beams),  Eqs.  (13)  and 
(14)  represent  the  two  light  bias  contributions:  Equation  (13)  defines  the  positive  definite  contribution 
on  the  line  that  contains  the  signal,  with  width  determined  by  the  finite  AOBD  aperture;  Eq.  (14) 
defines  the  raster  line  pattern,  with  each  reference  line  exhibiting  a  width  resulting  from  Fourier 
transformation  of  the  repeating  coarse  chirp  (in  Eq.  (C9),  summation  of  the  exponentials  produced  the 
array  of  delta  functions,  or  “comb”  filter). 
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The  detectable  signal  results  from  interference  of  the  signal  and  coarse  chirp  reference  legs. 
This  is  constructed  by  multiplying  Eqs.  (II)  and  (12)  and  their  complex  conjugates,  adding,  and 
integrating  over  the  duration  of  the  slow  chirp.  This  introduces  the  function  rect  (t/T2)  since  the 
integrand  contains  nonvanishing  dependences  on  time.  Making  the  substitutions  z  =  f/v,  t  =  y  /«, 
and  D  =  vT{  and  carrying  out  the  integration  yields 

(Signal  form  upon  detection) 

Es(x,  y)  =  (2 CT 2 / v)  exp{/ 4 7r/oy  /u]  exp{—i[2*p/T\  +  b]2 /4a\ 
exp[—iirbp/aT\)-  £  j  df  U(f)  sine  [*T\(avx  -  /)] 

p 

sinc[irT2(f  -  20y/u  -  p/T{)\  sinctTrTifavx  -  p/Tx)\,  (15) 

where  the  stationary  phases  of  Eq.  (CIO)  have  been  reintroduced  for  completeness.  The  three  sine 
functions  represent,  in  order:  (1)  a  coarse  spectral  estimate  generated  by  the  lens  Fourier  transforms 
of  the  signal,  yielding  a  peak  on  the  line  av:c  =  f,  (2)  a  fringe  modulation  envelope  in  the  y  direction, 
giving  a  peak  at  the  point  /  =  2/3 y  /u  +  p  /7'1 ;  and  (3)  a  grid  of  horizontal  reference  frequency  lines 
shifted  in  frequency  by  an  amount  proportional  to  y  by  the  slow  deflector.  The  signal  must  fall  on 
these  lines.  This  effect  is  adjusted  for  by  tilting  the  signal  beam  or  the  coarse  chip  beam  in  the  x 
direction  by  a  small  amount.  The  phase  exp{i4ir/0y  /«}  represents  the  usual  stationary  acoustic 
fringes  obtained  in  a  one-dimensional  time-integrating  process.  Thus,  Eq.  (15)  is  the  mathematical 
representation  of  the  2-D  output  signal  realized  at  the  detector  in  terms  of  the  Fourier  amplitude 
spectrum  U(f)  of  the  input. 
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